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Abstract 

A seminal milestone in lattice statistics is the exact solution of the 
enumeration of dimers on a simple-quartic net obtained by Fisher, 
Kasteleyn, and Temperley (FKT) in 1961. An outstanding related 
and yet unsolved problem is the enumeration of dimers on a net with 
vacant sites. Here we consider this vacant-site problem with a sin- 
gle vacancy occurring at certain specific sites on the boundary of a 
simple-quartic net. First, using a bijection between dimer and span- 
ning tree configurations due to Temperley, Kenyon, Propp, and Wil- 
son, we establish that the dimer generating function is independent of 
the location of the vacancy, and deduce a closed-form expression for 
the generating function. We next carry out finite-size analyses of this 
solution as well as that of the FKT solution. Our analyses lead to a 
logarithmic correction term in the large-size expansion for the vacancy 
problem with free boundary conditions. A concrete example exhibit- 
ing this difference is given. We also find the central charge c = —2 
in the language of conformal field theory for the vacancy problem, as 
versus the value c = 1 when there is no vacancy. 
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1 Introduction 



The problem of enumerating close-packed dimers on a finite simple-quartic 
net was solved by Temperley and Fisher fl], |] and by Kasteleyn in 
1961. An outstanding related but yet unsolved problem is the enumera- 
tion of dimers on a net with vacant sites 0]. Here, we consider this vacancy 
problem when a single vacant site occurs on the boundary. 

The difficulty associated with the vacancy problem is that, while the 
determinant whose square root yields the dimer generating function can be 
written down using the Kasteleyn formulation [Q], its evaluation is difficult. 
In 1974 Temperley reported an intriguing bijection relating close-packed 
dimer coverings with spanning tree configurations on two related lattices. 
This offers an alternate approach to the vacancy problem since spanning trees 
can be enumerated by standard means. The Temperley bijection has been of 
renewed recent interest and extended to graphs with certain weighted and/or 
directional edges Here, we use an extension of the Temperley bijection 
due to Kenyon, Propp, and Wilson to study the vacancy dimer problem. 

Our first result is that, for a simple-quartic net with free boundaries 
and one fixed vacant site located at certain specific sites on the boundary, 
the dimer generating function is independent of the position of the vacancy. 
The exact generating function for close-packed dimers on this net is then 
deduced from that of the spanning trees. In view of the connection with the 
conformal field theory and current interests in finite-size analyses for two- 
dimensional lattice models 0, ^ |T0| , we next carry out finite-size analyses for 
the weighted spanning tree solution as well as that of the FKT solution. It is 
found that a logarithmic correction term arises in the large-size expansion in 
the case of the vacancy problem with free boundaries, a term which is absent 
in the expansion of the FKT solution. A concrete example demonstrating 
this difference is given. We also find that the occurrence of a vacancy yields 
a new central charge c = —2 in the language of the conformal field theory. 

The organization of this paper is as follows: To make the paper self- 
contained, we restate and establish in section 2 the bijection due to Temper- 
ley between spanning tree and dimer configurations, as well as an extended 
version of the bijection due to Kenyon, Propp, and Wilson 0]. This extended 
Temperley bijection permits us to establish in section 3 the independence of 
the dimer generating function on the location of the vacancy when it occurs 
at specific boundary sites. The explicit expression of the generating function 
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for a simple-quartic net is then obtained. In section 4 we carry out finite-size 
analyses for weighted spanning trees and the FKT dimer solution. Specializ- 
ing the results to dimer enumerations, we find a logarithmic correction term 
which is unique to the vacancy problem with free boundaries. Discussions 
and a summary of our findings are given in section 5. 

2 The extended Temperley bijection 

For definiteness we consider a simple-quartic net (lattice) with free bound- 
aries, although much of the results of this section also hold for more general 
planar graphs @, |^ . 

First we restate the Temperley bijection. Starting from a Li x L2 simple- 
quartic net G with free boundaries, one constructs a dimer lattice by i) 
adding a new site at the midpoint of each edge of G, ii) inserting in each 
internal face of G a new site connected to the midpoints of the 4 edges of G 
surrounding it, and iii) removing one corner site of the resulting lattice and 
its incident edges on Gd. Thus, Go has a total of (2Li — 1)(2L2 — 1) — 1 sites 
consisting of the original L1L2 — 1 sites of G, which we call the odd sites, and 
the remaining (2Li — 1)(2L2 — 1) — L1L2 new sites, which we call the even 
sites. An example of this construction for Li = L2 = 3 is shown in Figs. 1(a) 
and 1(b). 

A spanning tree is a collection of connected edges of G which does not 
form closed circuits and covers all sites. Then we have the 

Temperley bijection: There exists a one-one correspondence between span- 
ning tree configurations on G and dimer configurations on Go- 
To see that the bijection holds, one observes that to each spanning tree 
configuration on G, one can construct a unique dimer configuration on Gd by 
first laying a dimer along each tree edge, starting from the edge(s) covering 
the corner site of Gd which has (have) been removed, and proceed along 
the spanning tree edges in an obvious fashion. After laying dimers along 
all tree edges, the remaining sites of G^) can then be covered by dimers in 
a unique way [Q. Conversely, starting from each dimer configuration on 
Gd, one constructs a unique tree configuration on G by drawing bonds (tree 
edges) along dimers originating from all odd sites. These bonds cannot form 
close circuits, since otherwise they would have enclosed an odd number of 
sites of Gd which is not permitted in close-packed dimer configurations. This 
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process leads to a unique tree configuration on G. This completes the proof. 
An example of the Temperley bijection is shown in Figs. 2(a) and 2(b). 

Kenyon, Propp, and Wilson p have shown that the Temperley bijection 
holds more generally for graphs with certain weighted and/or directed edges. 
For our purposes, however, we shall confine ourselves to the original Tem- 
perley bijection as stated in the above, with the step iii)replaced by one of 
removing an odd site on the boundary together with its incident edges. An 
example of constructing such a Gd is shown in Fig. 1(c). The proof of the 
bijection between tree configurations on G and dimer coverings on G/jgoes 
through as before, and we are led to the 

Extended Temperley bijection (Temperley-Kenyon-Propp- Wilson): There 
existsa one-one correspondence between spanning trees on G and dimer cov- 
erings on any Gd constructed from G by removing any boundary odd site and 
its incident edges in step Hi) of the construction described in the above. 

An example of such a bijection is shown in Figs. 2(a) and 2(c). 

Remark: The extended Temperley bijection does not hold for dimer lat- 
tices Gd containing an interior vacancy. In that case while each spanning tree 
can still be mapped into a unique dimer configuration as before, there exist 
dimer configurations which cannot be mapped into spanning trees. These 
are dimer coverings with no dimers laying on any of the 4 edges of G incident 
to the defect site. 

3 Dimer lattice with a vacant boundary site 

3.1 Dimer generating function 

The dimer generating function for is 

Z{Gn;x^,X2)= E x'l'xT (1) 

dimcr config. 

where the summation is taken over all dimer covering configurations, xi and 
X2 are, respectively, the weights of horizontal and vertical dimers., and rii and 
n2 are, respectively, the number of horizontal and vertical dimers. Clearly, 
we have 

Z{G£,]1,1) = the number of dimer configurations on Gd- (2) 
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Consider two different dimer lattices Go and G'j^ obtained from G as 
described in tlie above, namely, by removing different boundary odd sites. 
Then we have the following equivalence: 

Proposition 3.1.1: 

Z{Gd] xi, X2) = ^(G'^,; xi, X2) (3) 

for any Go and G'j^. 

Proof. The extended Temperley bijection dictates that there is a one- 
one correspondence between spanning tree configurations on G and dimer 
configurations on any Gn- It follows that there is a bijection between dimer 
coverings on Gd and G'^j, and that the summation in (|l|) on Gd can be 
considered as taken over all spanning tree configurations on G. 

For each spanning tree configuration T of G, the dimer weight in the 
summand in (0) consists of two factors, 

x'l^xr = Wo{T;x,,X2)We{T;x,,X2), (4) 

where Wo is the product of the weights of those dimers originating from odd 
sites, and We is the product of the weights of those dimers covering two 
even sites. For the two dimer coverings of Gd and G'j^ corresponding to the 
sameT, their We factors are the same by definition. Their Wo factors are also 
the same since, even though the respective dimer positions may be shifted, 
they lay along the same spanning tree edges hence carry the same weights. 
It follows that summations on the Ihs and rhs of are identical term by 
term, and the proposition is proved. • 

Remark: Proposition 3.1.1 holds more generally for arbitrary planar G 
and its related Gd- Since the overall dimer weight can always be factorized 
into the product WoWe as in (p, the proof of the proposition goes through 
as presented. 

We next consider the generating function of weighted spanning trees for 
the Li X L2 simple-quartic G. Assign weights xi and X2, respectively, to 
edges in the horizontal and vertical direction. Then, the weighted spanning 
tree generating function is 

TiG;x,,X2)=J2^T^T (5) 

T 
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where the summation is taken over all spanning tree configurations T on G 
and, as in (|l|), ni and n2 are the numbers of edges in the spanning tree in 
the respective directions. Particularly, we have 

T[G] 1, 1) = the number of spanning tree configurations on G. (6) 

From the extended Temperley bijection, it is clear that we can also write (|D 

as 

T{G;xi,X2) ='^Wo(T;xi,X2) (7) 

T 

where Wo is the factor in (^) for the dimer covering on any Go- It is seen 
from (|^ that if all dimers covering even sites of Go have the weight 1, namely 
We = 1, then the dimer generating function is simply T{G;xi,X2). 

More generally for a simple-quartic G of size Li x L2 and the related Gd, 
we have the equivalence: 

Proposition 3.1.2: 

V X2 xiJ 

Proof. From the construction of Gd we note that there are a total of 
L2{Li — 1) even sites located at midpoints of horizontal edges of G. We call 
these the H sites. Similarly, there are Li{L2 — 1) even sites of located at 
midpoints of vertical edges of G, which we call the V sites. 

The H and V sites can be covered by either horizontal or vertical dimers. 
Let Nfjh (Nhv) be the number of horizontal (vertical) dimers covering the H 
sites. Then we have 

NHh + NH. = L2{L,-l). (9) 

Likewise, we have 

Nvh + Nv, = L,{L2-l), (10) 

where Ny^ {Nyv) is the number of horizontal (vertical) dimers covering the 
V sites. In these notation, we can rewrite the spanning tree and dimer 
generating functions as 

T{G;x,y) = ^x^^'^y^- 

T 

Z{Gd; xi,X2) = 5]xl^«'^+^^'^a;2'^«''+^^^ (11) 

T 
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where we have used the one-one correspondence between spanning tree con- 
figurations T and dimer configurations. The Proposition 3.1.2 now follows 
after eliminating Nyh and Nhv in ( P^ using (||) and (|T0]). • 



3.2 Dimer enumerations 

For a simple-quartic G of size Li x L2 with free boundaries, the generating 
function for weighted spanning trees has been evaluated |Tl| and is given 
by 



T{G]Xi,X2) 



LiL 



1 L1-1L2-1 

n n 



2 xi[l — cos -;— ) + 2 X2 I 1 — cos — j 



2 m=0 n=0 I \ -^1 

(m,n)^ (0,0). (12) 

Now the dimer lattice Go is of size M x N with a boundary vacancy, where 

M = 2Li - 1, = 2L2 - 1. (13) 

In terms of the dimer lattice sizes M, A^, we thus have, after using Proposition 
3.1.2 and (O), 



Z{MxN-l}{xi,X2) 



X 



(Af-l)/2 (7V-l)/2 



M-1 
2 2 

X n n 

m=l n=l 



2 2 2 2 

4X1 COS — + 4X2 COS 



M+ 1 



Ar+ 1 



(14) 



Here, we must have MN = odd to admit dimer coverings. The subscript 
{M X — 1} in ([1^ reminds us that the enumeration is for an M x net 
with one boundary odd site removed. Note that the factor m = -n, = is 
excluded in the product in ([I^). This expression is to be compared with the 
enumeration of dimers on an M x A^ simple-quartic net without vacancies. 
For M and A^ both even, for example, the expression is M 



'{M,N} 



(Xi,X2) 



M/2 N/2 

n n 

m=l 71=1 



4x1 "^os^ 



mvr 



M+ 1 



+ 4x1 cos^ 



nvr 



A^ + 1 



(15) 



4 Finite-size analyses 

Finite-size expansions of physical quantities associated with two-dimensional 
lattice models have been of current interest both in physics 0, |, and 
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in mathematics [Q. For the dimer problem Kenyon has recently deduced 
very general results on the leading terms of the asymptotic expansion of the 
dimer enumeration for rectilinear lattices with free boundaries of any shape. 
Alternately, one can obtain expansions for regular lattices, in principle to all 
orders, by analyzing known exact expressions This is the approach we 
now use. 

4.1 Spanning tree generating function 



Consider first the generating function (|T2]) for the fully weighted spanning 
trees. For large Li and L2, we expect to have 

\nT{G]Xi,X2) = fhuik{xi,X2) + fc{xi,X2) (16) 



LiL 

where /bulk is the per-site bulk free energy and fc is the correction containing 
terms of the order of L^^, ^ and higher. Using (|1^), we find the bulk free 
energy 

/buik(xi,X2) = ^ lim -i^—\nT{G;xi,X2) 

Li,L2^co L1L2 

= — / de (i01n 2 xi(l - cos6') + 2 ^2(1 - COS0) 

71'^ Jo Jo J 

— / d(j) \n [J Xi + X2sm'^ (p + y/x2 sin (j)), (17) 



where the last line is obtained by carrying out the 6 integration. The com- 
putation of correction terms fc for products of the form of (|12|) is standard 
p|, H, [12| . Particularly, one has 



/buik(l,l) = -G (18) 

TT 

where G is the Catalan constant given by 

G = 1 - + 5-2 _ 7-2 + . . . = 0.915 965 594 ... . (19) 

To compute ([l6| ) we proceed as follows. Take out a factor xi from each 
of the L1L2 — 1 factors in(|l^) and split the product into 3 parts to take care 
of the exclusion of the m = n = factor. 
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We have 



T(G;xi,X2) = (L1L2) 'xfi^^-i (T0T1T2) 



where 



Li-l L2-I 



L2-I 



^0 = n n i^K^), ^1= n ^kq), n ^(o,^) 



m=l n=l 



m=l 



n=l 



with 



and 



mvr 



F(m, n) = 2 - cos + 2r - cos ^ j 
r = a;2/a;i. 



Using the identity |T3 



Af-l 

n 

n=l 



"nil 

2 cosh 2^-2 cos — 
N 



and its 6* — > hmit, 



Af-l 

n 

n=l 



2 - 2 COS ■ 



sinh2A^^ 
sinh2^ ' 



we find 



^ _ sinh(2Lift 
Jo - 11 



^^-^ sinh 26'„ 



T2 = r^^-^L2, 



with given by 



or, exphcitly, 



cosh 26n = 1 + T I 1 — cos — j . 



^n = F(^) = [cosh-i(l + 2a^)]/2 

= sinh~^ ttn 

= In (an + \/l + al^ , 



where 



ctn = VTsm 



nil 
2L. 



Substituting (^) into (|T^, we thus obtain 



(30) 



\ sinh 29 n 

The product in the denominator in (|30D can again be evaluated using (p^) as 



Yl sinh2 20„ 



n=l 



L2-I L2-I 

n (cosh20„ - 1) ■ n (cosh2^„ + 1) 

n=l n=l 

n Ki-^°^T-)1- n 



mr-- 

n=i -^2 ' - „=i 

2{L2-1) 



2 + r — r cos 



r 



/ sinh2L2a' 
\^ sinh 2a 



where a is given by 



cosh 2 a = l + 2r"^ 
or sinh a = 1/ -^/r, or exphcitly, 

a = In (Vt-^ + Vl + r-i^ 

Combining these results, we obtain from (^) the expression 



(31) 
(32) 
(33) 



sinh 2a; 



1/2 L2-1 



L2 sinh 2L2a^ 
Taking the logarithm, we obtain 



n [2sinh(2Li0„)]. (34) 



n=l 



Lo-l Lo-l 



lnT(G;a;i,X2) = (L1L2 - 1) Inxi + 2Li ^ ^„ + ^ ln(l-e-^^i' 

n=l n=l 

-L2 a - i In (1 - - - In L2 + - ln(sinh 2a). (35) 



For large Li and L2 with the ratio L1/L2 finite, the first two terms in 
([35| ) contribute to the bulk free energy /buik(a;i, 2^2) given in (|T7|) . To carry 
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out the summations in (P^, we use the Euler-MacLaurin summation formula 
given by 



N 



^ 1 fa+N5 I 



n=l 



+ ^[f\a + N5)-f{a)]+0{5'). 



(36) 



With a = 0, = La, 5 = 7r/2L2, and f{x) = F{x) defined in (|8D, one has 

L2—1 L2 

n=l n=l 

L 



[/buik(a;i,a:2) - lna;i - - ln(Vl + r + ^r) 



(37) 



2AL2 

For the second summation in (|3^) , we follow the manipulation in |P to write 
L2-1 



n=l 



-2nL\-K y/r I L'^ 



n=l 



+ 0(L, 



(38) 



for some < e < 2, with e ^ when L2 ^ 00. Putting the results together, 
we find the finite-size correction 



ri X Ci(xi,X2) , 02(0:1, X2) , C3(Xi,X2) , / 1 
fc{Xx,X2) = \ \ — h O 



LIL2 



(39) 



where 



Ci(Xi,X2) 
C2(xi, X2) 



In (Vt-^ + Vl + 



- In ( + vT+r 



-- lnL2 + ln2 + - In 5-^5 

2 4 xfxi 



cxj 

+ ^ In (1 - e-^^'^v^^i 



I2L2 



n=l 



(40) 



Particularly, for Xi = X2 = I, the expression /c(l, 1) given by (|39D reduces to 



the one given in M and |14 
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Despite its appearance, the expression for C3(xi, X2) is actually symmetric 
in {xi, Li} ^ {x2, L2}, a. fact can be seen from the identity 



ilnL,-llnx,-^^. 
2 4 I2L2 



m=l 

00 

n=l 

Introducing the Jacobi theta function 

00 

g) = 2gi/^ sin [] " " 2?'" cos 20 + g^"), 

n=l 



and the identity |15 



-In 



^'1(0, g) 



nl/3 



2gi/4 



n=l 

where "d'^ is the derivative of d\ with respect to 0, then we have also 
C3(xi, X2) = -i In L2 + ^ In ^^i^ + ^ In [4^;(0, g)" 

where 



-L\y/T-K j L-i 



The identity (^Tj) follows from the Jacobi transformation [|15 
(0, e""^) = (0, e-^l"), > 0. 



(41) 



(42) 



(43) 



(44) 
(45) 
(46) 



Now, the theta function 'd\ is finite for L1/L2 finite,^ it follows that the 
leading behavior of C3 is lnL2 (~ InLi). In conformal field theory ^ one 

needs to compute the limits 



1 lnT(Z2) 

— hm 

1 lnT(Z2) 

— hm 

1j2 il— >00 Lil 



/bulk + — + + ^(^2 )• 

L2 Ij2 



(47) 
(48) 



-'^This also follows from the fact that the series X^^^i ^^^(l ^ 
S^^i X]m=i w™"/"^ = ~ X]m=i "™/™(l ~ w™) converges for all < u < 1. 
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Using (PPP, we find 



Ai = A, = (49) 

This yields a central charge c = —2 upon taking xi = X2 = 1 (r = 1). 

4.2 Dimer enumerations 

We are now in a position to analyze the finite-size corrections of the two 
dimer enumerations (|T^ and (^3|). Although the expressions refer to two 
dimer lattices with different geometry, one for an M x — 1 lattice with a 
vacancy and MN = odd, and one for an M x lattice with MN = even, a 
comparison can still be meaningful if both expansions are expressed in terms 
of lattice sizes M and A^. 

(a) Close-packed dimers: For close-packed dimers on an M x A^ net with 
MN = even, we have carried out the analysis for the expression (^) along 
the lines outlined in the above, and obtained the result (which can also be 
extracted from discussions in [§,3) 

\nZ{MxN}ixi,X2) = (MN + l)/b,ik + Nci + MC2 + Cg + o(l), (50) 

where 

/buik(2;i,a;2) = ^/buik(a;i,X2), 



ci(xi,X2) = fbnik{xi, X2) - -In i^xi + ^ xf + xl 

C2(xi, X2) = /buik(a;i, X2) - ^ In (x2 + \lxl + x]^ 

C3(xi,X2) = ^ln2 - ^ In (^i + + ~\^^ (^^ + \l + 

+\ HA + 4) + ^ + E In (1 + e-(--)'^^-/^-) . 

(51) 

Expressions of Ci{xi,X2) and C2{xi,X2) reduce to those found by KenyonQ 
when Xl = X2 = 1- In the language of the conformal field theory [0], the term 
7tMx2/24Nxi in C3 yields the central charge c = 1 upon taken M = N and 
Xl = X2, the accepted value for dimer and Ising systems. 



13 



Again, the expression ( [^ ) for C3 is symmetric in {xi,M} ^ {x2,A^}, a 
fact can be seen from the identity 



24iVxi „tt 



-{2n-l)TTMx2/Nx-i 



TlNXi 



^ In (1 + e-(2™-l)'^Afa:'i/Ma;2 

(52) 

The series ln(l + u^^ ^) converges,^ so C3 does not diverge for large 

M,N. 

The expression for C3 can also be written as 



C3[Xl,X2) 



^ In 2 - ^ ln(xi + ^Jxl + xfj - ^ ln(x2 + + xl) 

+ ^ ln(x? + x^) + 1 In ^3(0, g) - i ln^9;(0, g), 
4 2 



where g = e '^^^a/A^aii g^j^^^ ^j-^g theta function ^2, is given by [15 

00 

^3(</>, g) = n (1 - ^'") (1 + V"-' COS 2</. + g 



,4n-2N 



(53) 



(54) 



ra=l 



Then, the identity (0) is a consequence of the Jacobi transformations (|T) 
and 

^?3(0, e-^'") = v^'^^^iiQ, e~"^), w > 0. (55) 



(^6^ Close-packed dimers with a boundary vacancy: For the simple-quartic 
net M X N — 1 with a boundary vacancy and MN = odd, one uses 



Li = (M + l)/2, L2 = (iV + l)/2, 
and expand ( pUD for large M and A^. After some algebra, we find 

In Z{MxN~l}{xi,X2) 

= (L1L2 - Li) Inxi + {L1L2 - L2) lnx2 + InTj G; — , — ) 

V X2 Xl 



(56) 



(L1L2 - Li) Inxi + (L1L2 - L2) lnx2 + L1L2 
{MN + l)/b,ik + Nci + Mc2 + c, + 0(1), 



Xl X2 



/bulk! — , — ) + fc[ 

\X2 Xl \X2 Xl 



(57) 



^This also follows from the fact that the series X]^^il^^(l + ^) 
X]m=i i—u)"^/m{l — u^"') converges for all < u < 1. 
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where Ci and C2 are given in (^Tj), and 



03(0:1, Xa) 



2 In 

2 2 



^ In (xi + /rfTxi) - ^ In (xs + 



j + ^ In ( 1 + 



2 12A^xi 



+ 5^ In (1 - e 



n=l 



'2n-KMxilNx\ 



Comparing C3 with C3 given by (^), we see that for large M, the deletion 
of a boundary site introduces a logarithmic correction term — In -y/iV in C3. 
Furthermore, upon taking M = N the term —7iMx2/12Nxi in C3 yields a 
central charge c = —2, which is different from that of the dimer system 
without vacancies. 

To verify the occurrence of a logarithmic term for the vacancy problem, 
we consider the following example. Consider two dimer nets of A^^ — 1 sites 
each, where iV > 3 is an odd integer so that the nets admit dimer coverings. 
While the two nets have different geometries, one a rectangular net of size 
(A^ + 1) X (A^ — 1) and one a square net of size N x N with one boundary 
odd site removed, they have the same area and perimeter. Any difference in 
the evaluationsof (0) and ( p7D would occur in C3 and Cg and higher order 
terms. Now from fPOl) and (157|) we obtain 



In Z 



{{jv+i)x{jv-i)}(l,l) = N^I^]+2Nci-\n{l + V2) + ^ln2 + ^ 



^In [l + e-(2"-i)- 



InZ 



{NxN-l} 



n=l 

r2 



+ 0(1) 



{N' + l) (^] +2Arci-ln(l + V2) + ^ln2 



TT 1 

12 ~ 2 



Defining the ratio 



R{N) 



In A^ + ^ In [1 - e" 

n=l 

^{(iV+l)x(iV-l)}(l, 1) 



2mT 



'{NxN-1} 



+ 0(1). (59) 



(60) 



and using (^91), we find the large A^ behavior 

R{N) CVn, N ^00, 



(61) 
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where 



R{N) 



i-G/w oo /i I p-(2n-l)7r' 



C = hm = TT , 



n=l 

0. 578 250 .. . . (62) 



As a numerical check, we have computed the value of R{N)/\/N for A^ = 3 
to 2251 using (|1^) and ([T5|). For A^ = 9, for example, one has 



Z|iox8}(l, 1) = 1 031 151 241 = (89) x (11 585 969) 
'{9x9-i}(l, 1) = 557 568 000 = 2^^ x 3^ x 5^ x (11)^ 

R{9)/V9 = 0.616 457... . (63) 



Results plotted in Fig. 3 confirm the large A^ limit of C given by ( p^ as well 
as the occurrence of the logarithmic correction in the vacancy problem. We 
remark that a similar result obtained by Kenyon 0] involving the occurrence 
of vacancy sites in the middle of a rectilinear net gives the ratio 

R{N) ~ N^/\ N^oo, (64) 

and thus a logarithmic correction — | In A^ in the finite-size expansion. 



5 Summary and discussions 

We have used the Temperley-Kenyon-Propp- Wilson bijection between span- 
ning trees on a lattice with free boundaries and dimer configurations on a 
related lattice with a boundary vacancy to establish the independence of the 
dimer generating function on the location of the vacancy. The equivalence is 
stated in Proposition (3). The generating function for close-packed dimers 
on a lattice with a single boundary vacancy is next computed, and compared 
with that of the known results for dimers without vacancies. It is found that 
the vacancy introduces a logarithmic correction in the large lattice expansion. 
A concrete example exhibiting this correction for an MN = odd net with a 
vacancy as compared to an MN = even net without vacancies is given. 

To ascertain whether the logarithmic correction is due to the defect of a 
vacancy, or due to the oddness of the net size, one needs to compare expan- 
sions for two nets (of the same even-even lattice), one with two boundary 
vacancies and one without vacancies. While this problem can be formulated 
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as the evaluation of the inverse of a matrix []TB[ in the Kasteleyn formula- 



tion, we argue that since the correction in question is that of the physical 
free energy of a dimer system, on physical ground one expects the correction 
to be additive for vacancies located sufficiently far apart. This would imply 
the logarithmic correction to be a "local" property due to the occurrence of 
vacancies. 

We have also found that in the language of the conformal field theory the 
central charge for the vacancy problem is c = —2 as compared to the value of 
c = 1 for the dimer solution without vacancies. Furthermore, the \fN ratio 
( [61| ) implies the existence of a boundary operator with scaling dimension 1/2, 
a value which does not appear in the standard Kac classification of operators 
at central charge — 2.|^ The extraction of the central charge should be viewed 
with caution, however, since the dimer systems do not exhibit critical points. 

We have also carried out finite-size analyses (details of which to be given 
elsewhere) for spanning trees on simple-quartic nets with other, including 
the toroidal, cylindrical, Mobius, and Klein bottle, boundary conditions. It 
is found that the logarithmic correction reported in this paper arises only in 
the case of free boundaries. This is consistent to the fact that the formulation 
of the extended Temperley bijection as presented in this paper is a property 
that is unique to graphs with free boundaries. 
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Figure captions 



Fig. 1. The construction of dimer lattices Go from a 3 x 3 spanning tree 
lattice G. Solid circles denote odd sites and open circles denote the odd site 
that has been removed in Gd- (a) A spanning tree lattice G. (b) A dimer 
lattice Gd constructed from G with one corner site removed, (c) A dimer 
lattice Go constructed from G with one odd site on the boundary removed. 

Fig. 2. The bijection between spanning trees on G and dimer configura- 
tions of Gd- (a) A spanning tree configuration on G. (b) The corresponding 
dimer configuration on the dimer lattice Gd oi Fig. 1(b). (c) The corre- 
sponding dimer configuration on the dimer lattice Gd oi Fig. 1(c). 

Fig. 3. The enumeration of R{N)/\fN for = 3 to 2251. The dashed 
line indicates the value C given by (^2|) in the large N limit. 
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